Quantifying Nonlocality Based on Local Hidden Variable Models 



o 
o 

(N 
> 

o 

as 



Oh 
I 



> 

ON 

(N 

G\ 
O 
ON 
O 



Dong-Ling Deng/ Jing-Ling Chen,^'* and Zi-Sui Zhou'^ 

' Theoretical Physics Division, Chern Institute of Mathematics, 
Nankai University, Tianjin 300071, People's Republic of China 
^College of Mathematics, Nankai University, Tianjin 300071, People's Republic of China 

(Dated: November 9, 2009) 

We introduce a fresh scheme based on the local hidden variable models to quantify nonlocality 
for arbitrarily high-dimensional quantum systems. Our scheme explores the minimal amount of 
white noise that must be added to the system in order to make the system local and realistic. 
Moreover, the scheme has a clear geometric significance and is numerically computable due to 
powerful computational and theoretical methods for the class of convex optimization problems known 
as semidefinite programs. 
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In a celebrated work in 1964, Bell showed surprisingly 
that different measurements on two separated parts of 
entangled quantum states would lead to stronger cor- 
relations, which are unexplainable by any local hidden 
variable (LHV) theory or shared randomness only [1]. 
This kind of stronger correlations is now well known as 
nonlocality. 

In the last decade, nonlocality has attracted much in- 
terest in both theoretical and experimental works [2] not 
only because of its close relations to the foundations of 
quantum mechanics, but also because of the vital role 
it plays in many quantum computation and informa- 
tion processes [3, 4]. Nevertheless, we are still far away 
from fully understanding nonlocality. Many fundamental 
questions still remain open, one of which is how to quan- 
tify nonlocality. Up to now, there are various schemes to 
quantify entanglement [5, 6], while few works for nonlo- 
cality quantification exist [7] . 

A nature consideration is to use the violations of Bell 
inequalities as the measure of nonlocality, based on which 
some nonlocality distillation protocols are proposed re- 
cently in the framework of generalized nonsignaling the- 
ories [8]. However, as these distillation protocols show 
exemplarily the advantages of this measure, they also 
drastically reveal its severe shortcomings: (i) For a spe- 
cific system, there are many tight Bell inequalities that 
are not equivalent to each other. Then which inequal- 
ity is the best one as the measure of nonlocality for this 
quantum system? (ii) Some quantum states do not vio- 
late a certain Bell inequality, but they do violate some 
other inequalities. What is worse, we even do not know 
whether these states violate Bell inequalities since we in 
principle need to find out all the Bell inequalities for the 
system to answer this question, (iii) Finding out all tight 
Bell inequalities for a quantum system is a NP hard prob- 
lem [9], so it is not pragmatic to use Bell inequality as 
nonlocality measure. Due to these disadvantages and the 
rapid theoretical and experimental developments in this 
realm, it becomes more and more pressing to introduce 
a more reasonable measure for nonlocality. 



In this Letter, based on the LHV models, we intro- 
duce a possible scheme with a clear geometric significance 
to quantify nonlocality for arbitrarily high-dimensional 
quantum systems. This scheme explores the minimal 
amount of white noise that must be added to the sys- 
tem so that the new system is local and realistic. In 
addition, by taking advantage of a brilliant work by Ter- 
hal et al. [10] and the established method for semidefinite 
programs, we show explicitly how to estimate efficiently 
this quantification by numerical method. 

To start with, we should specify some notations and 
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definitions. Let ^['^l = .^l 
note a rf-dimensional Hilbert space of n particle system. 
Here 3^^''^ , (fc = 1, • • • ,n), denotes the dfe-dimensional 
Hilbert space of the /cth subsystem and d — d\Y. ■ ■ ■ y. d^. 
Denoting the set of all quantum states on M'^'^^ that ad- 
mit LHV models by .^ihv, then the scheme to quantify 
nonlocality is as follows: For a specific quantum state 
p on the Hilbert space let = A-^ -I- (1 - A)p, 

(0 < A < 1), then the definition of the nonlocality is: 



M(p) — min A. 



(1) 



Here Id is a d-by-d identity matrix. Obviously, this defi- 
nition of nonlocality has the following properties: 

(i) M{p) is invariant under local unitary transforma- 
tions, i.e., M{p) =N{Ui (g) • • ■(S)UnpUl ®---®UX). 

(ii) M{p) = if the state p is separable. In fact, if p 
admits a HLV model, then it is obvious M{p) = 0. 

(iii) The J^{p) exists for any quantum state p since px 
is separable as long as A is large enough [11]. 

(iv) A/'(/ci) < 1. Actually, for any state p on J^I'^l, 
px is separable for A > 1 - 1 / ^/{d^^^T)id~Vj [12]. 
Thus we have a more accurate upper bound M{p) < 
l-l/V(d^-l)(d-l). 

(v) Af{p) > if p is an entangled pure state. This is 
obvious since any pure entangled state violates certain 
Bell inequahty [13]. 

Interestingly, this definition bears a geometric signif- 
icance as showed in the Fig. 1. It is obvious that the 
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FIG. 1: (Color online) A sketch to illustrate the geometric 
meaning of nonlocality. Here, the set of all quantum states is 
denoted by ^q; p* = N{p)^ + {l-N{p))p is the state with 
minimal A among all the states p\ that admit LHV models; 
pw ~ Id/d denotes the white noise state (random state). If 
we normalize the distance between nonlocal state p and the 
white noise state p^, then the nonlocality A/'(p) indicates the 
distance between p and p*. 

nonlocality M{p) indicates the distance between p and 
p* if we normalize the distance between nonlocal state 
p and the white noise state pw Moreover, from the ex- 
perimental point of view, M{p) also implies the minimal 
amount of white noise that must be added to the system 
in order to hide the nonlocal character of the state p. Its 
physical meaning is that it provides a lower bound of how 
much noise is needed. In other words, if the noise added 
to the system exceeds the bound, then the nonlocal cor- 
relations of this system will disappear. 

Now, we have the scheme to quantify nonlocality. How- 
ever, from the definition of JV{p), it is very difficult to 
compute J^{p) analytically since the analytical construc- 
tion of LHV models for some entangled states is itself ex- 
tremely difficult [14]. Fortunately, Terhal et al. have in- 
troduced a simple and efficient algorithmic approach for 
LHV-model construction for quantum states [10]. Their 
approach is based on the construction of a symmetric 
quasiextension of the quantum state. They showed ex- 
plicitly that a symmetric quasiextension might lead to a 
LHV model for the specific state, depending only on the 
number of local measurement settings for each observer. 
In the following, we shall show how to estimate numeri- 
cally A/'(p) by taking advantage of this method. We will 
focus on the two-qubit and two-qutrit systems. The gen- 
eralization to multipartite higher-dimensional systems is 
straightforward. 

For convenience and completeness, we first recapitu- 
late the main results of Terhal et al. [10]. Consider the 
following Bell-type scenario: Two distant parties, Alice 
and Bob, has a set of local measurements. Denote the 
number of measurements for Alice as Ala and let each 
measurement has Oa outcomes. Similarly, we denote the 
number of measurements for Bob as Mb and let each 



measurement has Ob outcomes. Moreover, let's define 

A 

where A: ®^ ^ ®A1 jg ^ permutation of spaces ^ 
in We say that p on ^ (g) ^ has a (A^^, A^^)- 

symmetric quasiextension when there exists a multi- 
partite entanglement witness Wp on Jff (g) Jlf^^" 
such that Tr^»A4a-i ^<sM^^iWp = p and Wp — Sa® 

Ssiy^p)- Terhal et al. showed that if p has a {M.a,.M.b)- 
symmetric quasiextension, then p definitely admits a 
LHV model when Alice and Bob have M.a and M.b arbi- 
trary measurements [10]. 

Now, the task of constructing LHV models is trans- 
formed to a new task of finding symmetric quasi- 
extension for a specific state. According to Ref. [10], 
the new task can be stated as a particular case of the 
class of convex optimization known as semidefinite pro- 
grams (SDP) [10], which correspond to the optimization 
of a linear function subject to a linear matrix inequality. 
Vandenberghe and Boyd [15] write the typical SDP as: 

. . . T 
minimize c x, 

subject to F(x) > 0, 

where c is a given vector of length and F{x) = 
Fq + J2i XiFi', Fi {i = 1, • • • are some fixed p-by-p 
Hermitian matrices. The vector x, also of length is 
the variable over which the minimization is performed. 
The SDP is said to be strictly feasible if there exists a 
vector X such that F(x) > is satisfied. The dual prob- 
lem corresponding to a SDP, also a SDP, reads: 

maximize — Tr[_foZ], 
subject to Tr[_FiZ] = Z > 0, 

where the Hermitian matrix Z is the variable over which 
the maximization is performed. Similarly, if there ex- 
ists a matrix Z > satisfying the trace constraints, the 
dual SDP is also said to be strictly feasible. A very im- 
portant relation between the primal and dual optimiza- 
tions was shown by Vandenberghe and Boyd [15]: // 
both the primal and dual forms of a SDP are strictly fea- 
sible, their optima are equal and achieved by some fea- 
sible pair (Tiopt, Zopt). The constructions of symmetric 
extensions correspond to the dual form of a SDP [10] 
with Fi = Sa Sb{K (Xi Ci = ri = Tr[Aip] and 
Fo = IdA '^Idg'^I- Here {A^} is the basis for the space of 
Hermitian matrices that operate on J^a ® ■ This ba- 
sis is orthogonal in the trace inner product Tr[AiAj] = Sij 
and Ao = IdA I VdAds- The corresponding v = 

{dAdB)"^ - 1 and = d^^d^". 

Based on the discussion above, in order to calculate 
Af[p], we only need to find the minimal A so that p\ ad- 
mits LHV model. Nevertheless, the above numerical con- 
struction of LHV model depends on the number of local 
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FIG. 2: (Color online) Numerical results of nonlocality, Bell- 
CHSH inequality violation and entanglement of formation for: 
(a) the two-qubit pure state IV") = sin6l|00) +cosS|ll); (b)tlie 
maximally entangled mixed states pmems- The nonlocality 
7V(p) is estimated by A/''^'^'(p) and the Bell-CHSH inequality 
violation is normalized by dividing a factor of 2\/2. 



measurement settings for each observer. Consequently, it 
might be more convenient to introduce a transitional def- 
inition of nonlocality also based on the number of local 
measurement settings for each observer: 



' (p) = min 



A, 



(3) 



where -5?/;^°'^''' is the set of all the states that ad- 
mit LHV models for (TWa, A^&)-settings. We term 
j^[M^.Mt](^p-^ the (AlQ,7W6)-nonlocality of state p. 
Mathematically, one has 



N{p) = 



Ma 



lim 



(4) 



From the relationship stated in Eq. (4), we can always 
use J\f^^'^'^''\p) to estimate M{p) as long as Aia and 
A4h are big enough. Actually, numerical results show 
that for many states, J\f^^'^'^''\p) vary very slightly 
with different {Ma, Mb)- So, we only need to compute 
j^[Ma,M,]^p'^ with small Ma and Mb- Consequently, the 
scheme to quantify nonlocality introduced above can be 
carried out numerically due to the efficient method to 
construct LHV models with (Ala, Alh)-settings. 

To illustrate explicitly how the method works, we 
give some examples here. We have implemented the 
mentioned corresponding semidefinite program using Se- 
DuMi [16]. The first example is about a set of two- 
qubit states: \t/j) — sin0|OO) -I- cos 61111). The nonlocality 
A/'dV')) with different 6 is showed in Fig. 2 (a). Note that 
here we use J\f^^'^\\ip)) to estimate A/'(|'0))- Also drawn 
are the entanglement of formation (EOF) [6] and Bell- 
CHSH [17] inequality violation of this set of states. From 
the figure, the maximal nonlocality occurs at 9 = 7r/4, 
i.e., the maximally entangled state |5') = -^(|00) -I- 111)), 
and the maximal nonlocality Af{\'^)) equals to 1/3. For 
two-qubit pure state, the nonlocality increases monotoni- 
cally as the EOF goes from to 1. Interestingly, based on 
the fact that Af{\'^)) = 1/3, we suspect that there might 



exist some two-qubit Bell inequalities with an improved 
visibility 2/3, which is stronger than the CHSH inequal- 
ity. However, we have not found such an inequality and 
the question proposed by Gisin is still open [18]. 

Another example involves the so called maximally en- 
tangled mixed states (MEMS), which have the maximum 
amount of entanglement for a given linear entropy [19]: 



//(7) 7/2\ 

1-2/(7) 



V7/2 /(7)/ 



(5) 



where 7(7) = 



The nonlocal- 



7/2, for 7>2/3 
1/3, for 7<l/3 
ity M {pmems), entanglement of formation £{pmems) and 
Bell-CHSH violation B{pmems) for this set of states with 
different 7 are plotted in Fig. 2 (b). From the fig- 
ure, the difference between our scheme and the one us- 
ing Bell inequality violation to quantify nonlocality is 
obvious: according to our quantification, the nonlocal- 
ity of Pmems arises when 7 > 0.35, while the Bell in- 
equality approach fails to quantify the nonlocality for 
0.35 < 7 < ^ since the state does not violate the CHSH 
inequality for this region. Moreover, it is worthwhile to 
note that £{pmems) > and N{pmems) — for the re- 
gion < 7 < 0.35, indicating that for this region, the 
state Pmems is entangled, but it admits LHV model. This 
provides us a fresh evidence from another aspect that en- 
tanglement and nonlocality are different. 

The third example concerns the two-qutrit system. 
Let us first briefiy introduce the CoUins-Gisin- Lindcn- 
Massar-Popescu (CGLMP) inequality [20] and the von 
Neumann entropy. For the two-qutrit system, the 
CGLMP inequality reduces to: 



[P(Ai -Bi)+P(Bi =yl2 + l) 
P(B2-Ai)]-[P(Ai=Bi-l 
P{A2 ^B2-l) + P{B2 = Ai - 



+ P{A2=B2) 

+ P(Bi = A2) 
1)] < 2. (6) 



Here we use the symbol = to denote equality modulus 
3 and denote the joint probability P{Ai = Bj + m) 
{i,j = 1,2) that the measurements Ai and Bj have out- 
comes that differ, modulo three, by m: P{Ai = Bj+ni) — 
J2l=Q^i^i = a,Bj = a + m). Another related con- 
cept, which aptly measures the entanglement of pure two- 
qutrit states, is the von Neumann entropy. For a pure 
state pab of two subsystems A and B, the von Neumann 
entropy of the reduced density operators is given by: 

y{pA) = y{pB) = -Tr[pAlnpA] = -TtIpbIupb], (7) 

where pA = ^^bIpab] and pB = T^^aIpab] are the re- 
duced density matrices. The state considered here is 
the generalized GHZ state of two-qutrit system: \(j>) = 
sin^sin/3100) -|-sin^cos/3111) -|-cos^l22). As for the two- 
qutrit case, we also draw the nonlocality, von Neumann 



4 














—A— S[i\iD> \ 












-^'V(|^>.,,) 















1 




' 1 1 

n/12 


1 

n/6 


1 

ji/4 


It/3 


1 1 1 
5JI/12 


jt/2 



Angle ^ 

FIG. 3: (Color online) Numerical results of nonlocality, Bell- 
CGLMP inequality violation and von Neumann entropy for 
the two-qutrit state — sin ^ sin/3|00) + sin^ cos /3|11) + 
cosf|22). The nonlocality A^d^)) is estimated by A/'[^'^'(|0)). 
The Bell-CGLMP inequality violation and the von Neumann 
entropy are normalized by dividing 4-\/2 and 2 In 3, respec- 
tively. In this figure, we have only plotted the curves with 
^ = 7r/6, ^ = 7r/3 and ^ = tt/2. 



LHV models for quantum states and the powerful nu- 
merical methods for SDP problems. To illustrate how the 
method works, we have provided some examples concern- 
ing two-qubit and two-qutrit systems. Our approach does 
not need any Bell inequality, thus it naturally circumvent 
the previously mentioned disadvantages of using Bell in- 
equality violation as nonlocality measure. This scheme 
sheds a new light on the quantitative understanding of 
quantum nonlocality. Interestingly, we have found that 
it may be very useful in quantum phase transitions [21]. 
It would also be interesting and significant to explore its 
possible applications in quantum information and com- 
putation science, which we shall investigate subsequently. 
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entropy and the Bell-CGLMP inequality violation with 
different ^ and f3 in Fig. 3. Note that the nonlocal- 
ity A/'(|(?!))) is estimated by J\f^^^^\\(j))). For the conve- 
nience of drawing, the Bell-CGLMP inequality violation 
and the von Neumann entropy of the reduced density 
operators are divided by 4a/2 and 2 In 3, respectively. 
In other words, B{\(j>)) = X^(|0))/4y2 and V(|0)) = 
r(|(/)))/(21n3). From the figure, M{\(t))), B{\(t))) and 
V(|0)) behave quite differently as ^ and /3 vary. In this 
case, the maximal nonlocality Af{\(t))) does not occur at 
the maximal entangled state. Actually, for the maximal 
two-qutrit entangled state |$) = ^(|00) + |11) + |22)), 

the nonlocality is A/'(|$)) = |, which is less than 1/2, the 
nonlocality of the state |0)^=7r/2./3=7r/4 = "Tf^l*^*^^ ^ \^^))- 

It is worthwhile to point out that for all the three ex- 
amples, the nonlocality Af{p) is estimated by Af^'^'^\p). 
Of course, one can use A/"!^"'^"! (p) [Ma, Mb > 2) as 
the approximation of JV{p) and the bigger A4a and Aib, 
the better the approximation. However, as Ma and Mb 
become bigger, it might need more time and EMS mem- 
ory to run the SeDuMi program on the computer. Actu- 
ally, our numerical results show that for the states under 
discussion, the variations of A/'''^"^'^''! (/o) with different 
Ma and Mb are very small. Thus, we believe Af^'^'^^p) 
is sufficient for a rough estimation ofAf{p). 

In summary, we introduced a reasonable scheme with a 
clear geometric significance to quantify nonlocality based 
on the LHV models. The scheme is numerically com- 
putable due to the systematic approach of constructing 
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